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Abstract. The detection and first identification of the 
deuterium Balmer emission lines, Da and D/3, in the core 
of the Orion Nebula is reported. These lines are very nar- 
row, have identical 11 kms -1 velocity shifts with respect 
to Ha and H/3, are probably excited by UV continuum flu- 
orescence from the Lyman (D i) lines and arise from the 
interface between the H II region and the molecular cloud. 
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1. Introduction 

Deuterium is believed to be entirely produced in the Big 
Ba ng an d then steadily destroyed by astration (Epstein et 
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al. 1976 ). Standard models predict a decrease of its abun- 
dance by a factor 2-3 in 15 Gyrs (e.g., Tosi et al. 1998| ). 
This picture is essentially constrained by deuterium abun- 
dance determinations at ~ 15 Gyrs (primordial intergalac- 
tic clouds), 4.5 Gyrs (protosolar) and 0.0 Gyrs (interstel- 
lar medium). Although the evolution of the deuterium 
abundance seems to be qualitatively understood, the mea- 
surements show some dispersion. Thus, absorption in the 
Lyman series provides interstellar deuterium abundance 
(D/H) /SM ~ 1.5 x 10 -5 (Linsky |1998| ), but with fluctu- 
ations that may well be real (Vidal-Madjar et al. 1998). 
These dispersions led to the development of non-standard 
models in which, for example, deuterium may either de- 
crease by more than a factor 4 in 15 Gyrs (e.g., Vangioni- 
Flam et al. 1994) or be created/destroyed by new mecha- 
nisms [e.g., Lemoine et al. ( 1999] ) for a review]. 

A detailed appraisal of the evolution of deuterium is 
crucial for cosmology and galactic chemical evolution. The 
most reliable estimate of (D/H) /SM to date is based on 
far-UV observation from space (Copernicus, IMAPS, HST 
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or FUSE) of the Lyman lines of D and H in absorption. 
These lines are also observed in the optical and near-UV 
to obtain D/H in high redshift quasar absorbers. Other 
D/H determinations include in situ measurements in the 
Solar System (e.g., Mahaffy et al. 1998), observations of 
molecules such as HD or DCN (e.g., Bcrtoldi et al. 1999) 
and observations of D I 92 cm (e.g., Chengalur et al. 1997). 

New methods to determine D/H are of interest. One 
possibility is ground-based observation of the deuterium 
lines. The isotope shift of the deuterium Balmer lines with 
respect to the hydrogen Balmer lines is —81.6 kms -1 . 
These D I lines have never been identified before. Atte mpts 
to detect Da in absorption in the Sun (Be ckers 1975 ) and 
early-type stars (e.g., Vidal-Madjar et al. 1988 ) were un- 
successful (D is destroyed in stars). Traub et al. ( 1974 ) ob- 
served Ha in the Orion Nebula using three-etalon Fabry- 
Perot spectrometers and reported D/H upper limits. 

Here we report on spectra of Orion, secured at the 
Canada- France-Hawaii Telescope (CFHT). Emission lines 
detected in the blue wings of Ha: and H/3 are identified 
with Da and D/3. A preliminary account was presented 
by Hebrard et al. (1999). Observations are described in 
Sect. 0, the identification and the origin of the lines in 
Sect. |3[and ^ and the excitation mechanism in Sect. |5|. 

2. Observations, data reduction and results 

Observations of the Orion Nebula (M 42, NGC 1976) were 
conducted at the 3.6m CFHT, using the Echelle spectro- 
graph Gecko at the Coude focus with a slit length of ~ 40". 
The Ha and H/3 spectral ranges were observed in Octo- 
ber 1997 and September 1999 respectively. For Ha, the 
entrance slit was 1.2mm wide (3.5" on the sky), providing 
a resolution R = A/AA ~ 40 000 (~ 7.5 kms -1 ); the de- 
tector was the 2048 x 2048 "Loral 5" thin CCD and the 
spectral range was 6544A - 6576A. For H/3, the slit was 
0.8mm wide (2.3") leading to R ~ 50000 (~ 6 kms -1 ); 
the detector was the 2048 x 4500 "EEV2" thin CCD and 
the spectral range was 4832A - 4885A. 

The slit was centred 2.5' South of 9 1 Ori C (HD 37022), 
the brightest star of the Trapezium. The slit orientation 
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was slowly rotating during the exposures (Coude focus). 
Totals of 4.5 and 1.5 hours were devoted to Ha and H/3 
respectively, divided in 30 - 45 min sub-exposures. Small 
rotations of the grating were applied between Ha sub- 
exposures in order to disclose ghosts that may depend on 
grating setting. Ha was also observed at higher resolution 
(R ~ 80 000) in the same area for 20 min. Finally, H/3 was 
observed 2.5' North and 20" South of 9 1 Ori C with shorter 
exposures and R ~ 50 000. Bias, flats and Thorium-Neon 
lamp calibration exposures were secured regularly during 
the observations for each instrument configuration. 

The spectra were reduced using MIDAS software. The 
steps of the data reduction were as follows: (1) bias sub- 
traction; (2) flat division; (3) bad pixel and cosmic clean- 
ing; (4) summing the rows to transform the 2D-spectra 
into ID-spectra; (5) wavelength calibration; (6) shift to 
the heliocentric frame; (7) alignment of the different sub- 
exposures; and (8) suming up of the sub-exposures. After 
shifting to the heliocentric frame, both Ha and H/3 were 
fitted by a Gaussian on each sub-exposure. The standard 
deviation of the Gaussian peaks was less than 1 kms -1 . 
Sub-exposures were shifted to the average peak before 
summation in order to preserve the spectral resolution. 
An interfering signal, instrumental in origin, appeared in 
the H/3 spectra, producing small oscillations in the dis- 
persion direction, which slightly increased the noise level. 
Wavelengths are determined to better than 1.5 kms -1 and 

I. kms -1 in the Ha and H/3 final spectra respectively. 

Two weak emission lines are obvious in the blue wings 
of Ha and H/3 (Fig |l|). Anticipating the conclusion of 
Sect. ||, the weak lines are already identified with Da and 
D/3 in Table [j], where results of Gaussian fits to Ha, H/3, 
Da and D/3 are given. The full widths at half maximum 
(FWHM) of the deuterium lines are much smaller than 
those of the hydrogen lines. Relative fluxes (last row of 
Table [j]) are based on the theoretical Ha/H/3 ratio, thus 
implicitly correcting for reddening. For an H + -weighted 
electron temperature (0.85±0.10) x 10 4 K arid electron den- 
sity ~ 5xl0 3 cm -3 (e.g., Esteban et al. 1998 ), the Case B 
recombin ation ratio /(Ha)//(H/3) is 2.91±0.03 (Storey & 
Hummer 1995 ). Departure of this ratio from Case B is ex- 
pected to be much less than 1% in this thick n ebula for 
any reasonable dust content (Hummer & Storey 1992 ). 

Da and D/3 are seen all along the 40" slit. D/3 is present 
at all three positions. The velocity shifts between D/3 and 
H/3 at 2.5' N, 20" S and 2.5' S of 9 1 Ori C are respectively 

II. 8, 9.1 and 10.0 kms" 1 and the D/3 fluxes 4.2 ± 1.1, 
2.3 ± 0.6 and 5.7 ± 1.1 (H/3= 10 000). 



3. Identification of Da and D/3 

According to Table |l], the shift of both weak lines with re- 
spect to the hydrogen lines is —71 kms" 1 whereas the iso- 
topic shift of deuterium is —81.6 kms" 1 . This significant 
difference forces us to consider alternative possibilities. 
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Fig. 1. Vicinity of Ha and H/3 in Orion, showing Da 
and D/3 in emission. The x-axis is in kms -1 relative 
to the rest wavelength of either Ha or H/3. The verti- 
cal scale corresponds to peak fluxes 7250 and 2500 for 
Ha and H/3 respectively. Then /(Ha)//(H/3) = 2.91 and 
/(Da) //(D/3) ~ 1.10. 



• Spectral artifact? No feature equivalent to these lines 
is present in either the red wings of Ha and H/3 or the 
wings of [N n] A6548.05A, observed simultaneously. Small 
rotations of the grating applied between sub-exposures re- 
sulted in no change in profile, position and intensity of the 
lines. Finally no such lines were detected in bright plane- 
tar y neb ulae we observed in September 1999 (Hebrard et 
al. 2000 ), using the same instrument. This also excludes 
possible sky-line emission. In fact, no sky-lines have been 
reported at these wavelengths. Thus, these lines are real 
features, specific to the Orion Nebula. 

• Unidentified process or element? Attempts to find 
other identifications were unsuccessful. These lines can- 
not be scattered stellar emission [for example, Ha from 



Ori C is variable (Stahl et al. 1996)], as lines from 



hot stars are broad. For the same reason, they cannot be 
Raman features. The wavelengths do not correspond to 
any known quasi-molecular line. The fact both lines have 
identical velocity shifts with respect to H I practically re- 
stricts the possibilities to H I and D I emission (no He II is 
detected in the Orion Nebula). 

• High-velocity hydrogen emitting structure? Traub 



et al. (1974) reported the detection of a line in the blue 



wing of Ha. This line may correspond to ours, although 
it was interpreted by these authors as high- velocity Hi 
emission, noting the existence of a similar component in 
[Om] (Dopita et al. |1973[). Indeed, our R = 80 000 spec- 



trum shows a blue component in [Nil] but with velocity 
shift only ~ —22 kms" 1 . More importantly, any compo- 
nent arising from the H II region should have a minimum 
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Table 1. Gaussian fitted line profiles 



Line identification 


Ha 


Da 


H/3 


D/3 


Rest wavelength (A) 


6562.796 


6561.010 


4861.325 


4860.003 


Observed wavelength (A) 


6563.12 ±0.03 


6561.59 ±0.03 


4861.60 ± 0.02 


4860.44 ± 0.02 


v e (kms" 1 )" 


14.8 ± 1.5 


26.5 ± 1.5 


17.0 ± 1.0 


27.0 ± 1.5 


FWHM (kms -1 ) 6 


32.0 ±0.5 


8.6 ± 1.0 


32.1 ± 0.5 


8.1 ± 1.5 


Relative flux c 


29 100 ± 300 


6.3 ±0.6 


10 000 ±200 


5.7 ± 1.1 



a Heliocentric velocity. 

b Full width at half maximum corrected for instrumental width (original FWHM: 32.9, 11.4, 32.7, 10.1 kms -1 respectively). 
c Using the theoretical 7(Ha)/l(H/3) (see text), then J(Da)//(D/3) = 1.10 ± 0.22. 



width corresponding to thermal broadening. The thermal 
FWHM for hydrogen at 8500K is 20 kms -1 , much larger 
than 8 kms -1 (Table |l|). The -22 kms -1 ionized hydro- 
gen emission should be lost in the Ha and H/3 wings. 

Nonetheless, the fluorescence mechanism proposed be- 
low to explain the D I line excitation (Sect. ^J) may a priori 
apply to H I as well. One cannot formally exclude H I flu- 
orescence emission from a neutral, cold (thermal FWHM 
is ~ 7 kms -1 at 10 3 K), high-velocity (-74 kms -1 LSR), 
low velocity dispersion (<< 10 kms -1 ) layer keeping the 
same kinematical properties over many arc minutes. How- 
ever this is very unlikely as this hypothetical structure 
should in addition have a small column density (no other 
fluorescent line is seen at that velocity) and lie sufficiently 
close to the Trapezium stars (fluorescence varies as the in- 
verse square of the distance to the continuum source) . The 
survival of a neutral thin shell against photoionization (no 
low-ionization material is interposed; see Sect. 4]) is also 
in question. As a matter of fact, Cowie et al. (1979) de- 
tected several components of high-velocity gas in absorp- 
tion against l Ori (a star located within half a degree of the 
region we observed), notably a component at —68 kms -1 , 
close to —74 kms -1 . According to these authors, this com- 
ponent corresponds to a very old highly ionized supernova 
remnant situated over 100 pc from the Trapezium, at any 
rate too far away to yield fluorescence. 

It can therefore be safely concluded that these lines 
are Da and D/3. Kinematics (Sect. ^) brings out one more 
fundamental piece of evidence. 

4. Origin of the lines 

In the Orion Nebula, the Hn region is essentially matter 
bounded toward the observer and radiation bounded in 
the opposite direction (e.g., Rubin et al. 1991). The nar- 
rowness of the D I lines implies that they must originate 
in a cold, localised region along the line of sight, that is 
behind the H + region, in the "Photon Dominated Region" 
(PDR) where deuterium is in atomic form. 

This is borne out by available information on ve- 
locities. The heliocentric velocity of the Di lines is ~ 
27 kms -1 compared to ^ 16kms -1 for Hi. In a b lue spec- 
trum of the Trapezium region (Kaler et al. 1965|) , the Sill 
lines 3856±63A, probably produced by fluorescence in the 



PDR, appear shifted by ±11 kms -1 relative to the neigh- 
bouring H I Balmer lines, a shift similar to the one found 
for Da and D/3. Over a regio n close to the one we observed, 
Esteban & Peimbert ( 1999 ) measured heliocentric veloci- 
ties 6.4±1.4, 14.0±2.0, 24.6±2.2 and 26.8±1.4 kms -1 for 
Ar ++ , H + , O and N° respectively, tracing the free ex- 
pansion of the H II region, moving away from the molec- 
ular cloud. Over the same region, Hanel (1987) found 
20-23 kms -1 for [Nil] [in agreement with our measure- 
ment u([Nll]) ~ 21 kms -1 ] and 23.5-28 kms -1 for [Sn], 
thus encompassing the velocity we found for Da and D/3. 
From millimeter and submillimeter observations, Hoger- 
heijde et al. (1995) found velocities ~ 28 kms -1 for dif- 



ferent molecules. Observations thus clearly imply that the 
D I lines could arise from the boundary of the H n region. 

5. Fluorescent excitation of Da and D/3 

The narrowness of Da and D/3 allows the possibility to 
be excluded that the lines are excited by recombination 
in an ionized gas (Sect. ||). The H° column density of 
the PDR is - 10 22 cm -2 (Tielens k Hollenbach |l985[ ), so 



the D° column density is ~ 10 17 cm 2 [assuming a typical 



(D/H) 



ISM- 



10 ] and the optical thickness in, e.g., Lv/3d 



is over 100. Since the Di emission is confined to a layer 
coincident in velocity with that of the PDR (Sect. ||) and 
since the dust opacity there is orders of magnitude less 
than the Ly/3n opacity, fluorescence from the Lyman lines 
is a viable process to produce the deuterium Balmer lines. 
The UV continuum is dominated by 9 1 Ori C, whose ef- 



fective temperature is close to 4x 10 4 K (Rubin et al. 1991). 
Let us assume that both the ionization of the H II region 
and the deuterium fluorescence are due to a 4 x 10 4 K 
black body and that half the ionizing photons escape from 
the Hn region in the matter bounded directions (Ru- 
bin et al. 1991). If each photon impinging on the PDR 
at the Ly/3n wavelengths ultimately produces a Da pho- 
ton by scattering on D°, and only Lv/3d photons lead to 
Da excitation (neglecting cascades), then the flux ratio 
I(Dct)/I(Ha) is about 1.5 x 10 -4 x (Au/5kms -1 ), where 
Av is the full velocity width of the zone where Da is 
effecti vely excited. According to Tielens and Hollenbach 
( 1985 ) , the turbulent pressure in the PDR corresponds to 
Av ~ 5 kms -1 and according to Table 1, Av is proba- 
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bly less than 8 krns" 1 . The rather good agreement of this 
very coarse estimate with the observed value 2.2 x ICR 4 
(Table 1) is fortuitous. This estimate may be wrong in 
different ways. Many LyD lines can a priori absorb pri- 
mary photons and feed Da by cascades, then leading to 
an overestimation. Conversely, part of the LyD photons 
are absorbed by dust and/or reflected back to the Hn re- 
gion. Also, our estimate is global in character and our par- 
ticular line of sight may not intercept identical fractions 
of the Hn region and the PDR. Most importantly, the 
stellar continuum may be depleted in the vicinity of the 
Lyn lines, particularly for the first members of the series. 
Nonetheless, since these different effects tend to partially 
compensate one another, the above agreement indicates 
that the assumption of UV continuum fluorescence leads 
to the correct order of magnitude for the Da flux. 

Unlike Da/Ha, the flux ratio Da/D/3 is little sen- 
sitive to aspect effects. Assuming a ratio of visual ex- 
tinction to column density of hydrogen nuclei Ay / Nn = 
5xl0 -22 mag cm 2 and Apuv/Ay = 5 (Tielens and Hol- 



can provide independent information on the primary con- 
tinuum flux and on the competition of line scattering with 
dust absorption for photons. O I fluorescence lines have 
been detected long ago in Orion and the excitation process 



lenbach 1985), the ratio of Ly/3o opacity to dust opacity is 



~ 240. Only for large principal quantum numbers should 
dust absorption decrease D I fluorescence (the photoexci- 
tation cross section goes roughly as n~ 3 ). Since the stellar 
continuum should be about flat over the Lyman line range 
(except possibly in the vicinity of strong lines), some in- 
sight into the excitation process can be gained by assum- 
ing that all Lyo(n) lines, with principal quantum number 
n up to some given no, convert identical numbers of UV 
photons by fluorescence and that no fluorescence occurs 
for n > hq. Then, working out the cascades, the theoreti- 
cal Da/D/3 flux ratio is about 1.28, 1.41 and 1.51 for no = 
4, 5 and 6 respectively and tends to level off for larger 
no's. Only for no = 4 is this simple description compati- 
ble with the observed value 1.10±0.22 (Table 1). Since one 
would expect that a relatively large number of Lyn lines 
should contribute to the excitation, the suggestion is that 
one of the above assumptions was oversimplified or some 
significant process has been overlooked. For example, the 
stellar continuum is probably depleted in the wings of the 
Lyn lines and the vicinity of Ly/3jj is likely to be most 
affected, then selectively reducing the Da emission. Ob- 
serving higher deuterium Balmer lines is essential before 
attempting any detailed modeling. 

Note that the Da/D/3 of Table @ was obtained assum- 
ing that the reddening correction was the same for the H I 
and D I emitting zones. If extinction internal to the nebula 
is significant, then the actual (de-reddened) Da/D/3 ratio 
will be even smaller since the PDR, where the D I lines 
come from, is more deeply embedded than the H II region. 

Dust absorption will dominate LyD fluorescence for 
sufficiently large n, the deuterium Balmer decrement then 
changing from very flat to very steep. This break can lead 
to a D/H value inasmuch as the dust opacity per hydrogen 
nucleus is known. On the other hand, fluorescence lines 
from species co-extensive with D° including O I and Si II 



was established by Grandi (1975). Comparing Di and Oi 
lines may lead to a D/O abundance ratio. Detailed obser- 
vations and proper modeling of many deuterium Balmer 
lines and other fluorescence lines appear as a potentially 
accurate means to determine D/H in Hn regions. 

6. Conclusions 

Deuterium is identified for the first time in a nebula from 
optical spectroscopy. The excitation mechanism of the ob- 
served lines, Da and D/3, is continuum fluorescence from 
Lyo lines in the PDR. Considering the saturation of the 
first Lyo lines and the possible influence of the neighbour- 
ing Lyn lines, observing the full deuterium Balmer series 
is essential to obtain a D/H value from optical data and 
appears feasible, at least in Orion. 

An optical determination of D/H in H n regions would 
allow to check existing (D/H) JSM values and obtain 
D/H in low-metallicity extragalactic Hn regions, where 
the deuterium abundance should be close to its primor- 
dial value. 

The large photoexcitation cross section of the first Ly- 
man lines makes deuterium Balmer fluorescence a sensi- 
tive way to detect deuterium in nebulae, leading for exam- 
ple to stringent upper limits to D /H in planetary nebulae 
(Hebrard et al. 200C), where deuterium is depleted. 
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